Location of crossings in the Floquet spectrum of a driven two-level system. 
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Calculation of the Floquet quasi-energies of a system driven by a time-periodic field is an efficient way 
to understand its dynamics. In particular, the phenomenon of dynamical localization can be related to the 
presence of close approaches between quasi-energies (either crossings or avoided crossings). We consider here 
a driven two-level system and study how the locations of crossings in the quasi-energy spectrum alter as the 
field parameters are changed. A perturbational scheme provides a direct connection between the form of the 
driving field and the quasi-energies which is exact in the limit of high frequencies. We firstly obtain relations for 
the quasi-energies for some common types of applied field in the high-frequency limit, and then show how the 
locations of the crossings drift as the frequency is reduced. We find a simple empirical formula which describes 
this drift extremely well in general, and which we conjecture is exact for the specific case of square-wave 
driving. 
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I. INTRODUCTION 



The two-level system is a simple model which has been ap- 
plied to a great variety of physical problems. One application 
of growing importance is that of quantum computing 1 , since 
any quantum two-level system has the potential to act as a 
quantum bit. For this reason the coherent control of quantum 
states in these systems has recently become the focus of in- 
tense investigation 2 - 3 . A concrete example of such a two-level 
system is provided by a particle tunneling between two poten- 
tial wells, which can be experimentally realized by confining 
an electron to a pair of coupled quantum dots 4 . One method 
of controlling such a system, without destroying its coherence, 
is to apply oscillatory electric fields 5 . Such fields are able to 
produce the phenomenon known as coherent destruction of 
tunneling (CDT), in which the tunneling of the particle is sup- 
pressed when the parameters of the field are tuned to various 
"magic" values. As the applied field is time-periodic, Floquet 
analysis 6 has been applied to explain this non-intuitive result, 
and it has been shown 7 ' 8 that CDT is closely related to the 
presence of crossings or avoided crossings in the spectrum of 
Floquet quasi-energies. 

The driving field most frequently considered is of sinu- 
soidal form, and studies using CDT as a means of quantum- 
control have generally concentrated on varying either the 
envelope 9 or frequency 10 of a sinusoidal signal. In this work, 
however, we instead consider the effect of altering the signal's 
waveform. By using a perturbational method we first show 
how the waveform can be directly related to the quasi-energy 
spectrum, and give analytic results for sinusoidal, square- 
wave and triangular waveforms. These results are precise 
in the limit of high frequency. As the frequency is reduced, 
however, the locations of the crossings drift away from these 
values. This effect is extremely difficult to treat analytically, 
and such efforts 111213 produce complicated results which are 
difficult to interpret. Empirically, however, we find a simple 
formula which describes the drifting with good accuracy for 
many waveforms, and appears to be exact for the case of the 
square-wave. We thus provide a means for predicting the lo- 
cations of quasi-energy crossings for a given driving field in 
both high and low frequency regimes. 



II. METHOD 



A. Model Hamiltonian 



We consider a charged particle confined to a double 
quantum-dot system, described by the Hamiltonian: 



H = t{c\c R + h.c.) + (E L (t)n L + E R (t)n R ) 



(1) 



where the subscript L/R denotes the left/right quantum dot, 
ct and c . are creation and annihilation operators for a particle 

in dot /', and nj = c\cj is the usual number operator. The 
tunneling between the two dots is described by the hopping 
parameter t, and Ej(f) is the electrical potential of the external 
driving field. Clearly only the potential difference between 
the two dots is physically of importance, and so we can use 
the symmetric parameterization: 



E 

El = ^f(t), 



(2) 



where E is the potential of the driving field and fit) is a T- 
periodic function describing its waveform. The Hamiltonian 
has been written using a basis of localized states, but it 
may be easily transformed to the standard two-level form via 
a SU(2) rotation, yielding the result: 



A E 
H = -o- z + -f(t)o- x , 



(3) 



where <x, are the standard Pauli matrices. In this represen- 
tation the basis states used are extended states, formed by 
symmetric and anti-symmetric combinations of the localized 
states. In the absence of a driving field [E = 0) it is clear that 
the two eigenstates of this Hamiltonian consist of a symmetric 
ground state, and an excited anti-symmetric state. The split- 
ting between these two levels is given by A, which is related 
to the inter-dot tunneling via A = 21. 



B. Floquet theory 

As the function f(t) is periodic in time, the Floquet the- 
orem may be used to write solutions of the time-dependent 
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Schrodinger equation as i//(t) — exp[~2e ; ?]0;(f), where <pj(t) 
is a function with the same periodicity as /(f) and is called a 
Floquet state, and £j is termed the quasi-energy. Although the 
Floquet states have an explicit time-dependence, their period- 
icity means that the dynamics of the system on time-scales 
much larger than the period of the driving field is effectively 
given only by the quasi-energies. In particular, if the two 
quasi-energies approach degeneracy, the dynamics of the sys- 
tem on this time-scale will appear to be frozen, producing the 
effect of CDT. Consequently, determining the quasi-energies 
provides a simple and direct way of studying the long time- 
scale behavior of the system, and indicates whether CDT can 
occur 14 . In this work we restrict our attention to driving func- 
tions which possess the symmetry fit) = -fit + T/2). Impos- 
ing this restriction means that the Hamiltonian Q is invariant 
under the generalized parity operation x — > -x, t — > t + T/2, 
and as a consequence the two Floquet states will also possess 
this symmetry, one being even and the other being odd. The 
von Neumann-Wigner theorem 15 thus allows the two quasi- 
energies to cross as an external parameter, such as the field 
strength, is varied. Breaking this symmetry by choosing an 
alternative form for the driving field would mean that the 
quasi-energies would be forbidden to cross, and thus close 
approaches between the quasi-energies could only consist of 
avoided crossings. 

The Floquet states and their quasi-energies may be conve- 
niently obtained from the eigenvalue equation: 



d 

Hit)-l- 



<jit) = erfjii). 



(4) 



To obtain approximate solutions to this equation we follow 
a perturbation scheme introduced originally by Holthausi& to 
treat both the two-level system and driven superlattices, and 
which was generalized recently to also include the effects of 
inter-particle interactions 17 . In this approach the Hamiltonian 
{3 is divided into two parts: H, which contains the tunnel- 
ing terms, and Hj which holds the electric field terms. We 
then find the eigensystem of the operator 'Hiit) — Hi — id/dt 
by working in an extended Hilbert space of time-periodic 
functions 18 , and apply the tunneling Hamiltonian as a pertur- 
bation. A consequence of dividing the Hamiltonian in this 
way is that the perturbation theory works well in the high- 
frequency limit co » t , but breaks down in the opposite limit 
when the tunneling provides the dominant energy-scale of the 
problem^. 

For the Hamiltonian given in EqQ] the problem of finding 
the eigensystem of Iriiit) simply requires the solution of two 
uncoupled differential equations: 



E d\ 
E d\ 



(5) 
(6) 



These can be integrated immediately, giving the solutions: 

4> ± it) = exp[+/£F(f)/2]exp[/e ± f], (7) 

where F(t) = f fit')dt'. (8) 
Jo 



The periodicity of the Floquet states clearly requires that e ± = 
mod co. Without loss of generality we can restrict the quasi- 
energies to lie in the "first Brillouin zone" i-co/2 < e < co/2), 
and thus to lowest order in the perturbation theory they are 
degenerate and zero. Standard degenerate perturbation the- 
ory can now be used to evaluate the first-order correction to 
the quasi-energies, requiring only that we work in the the ex- 
tended Hilbert space of r -periodic functions by defining an 
appropriate scalar product: 



«0ml0n»; 



{4> m it')\(t> n it'))dt' 



(9) 



where (-|-) is the usual scalar product for the spatial component 
of the wavefunctions, and (-|-)r denotes the integration over 
the compact time coordinate. 

As the tunneling component of the Hamiltonian H, is acting 
as the perturbation, the first-order approximation to the quasi- 
energies is given by the eigenvalues of the perturbing matrix: 



«ff«»7 





n<fih 



t(4> 2 -h 
o 



(10) 



Comparing this expression with the original tunneling Hamil- 
tonian Q reveals that the action of the applied field is to renor- 
malize the tunneling terms by the factors (4>1)t- As cf> + is the 
complex conjugate of the quasi-energies take the simple 
form: 



4i 



where 



>7 



{4> + )t = - exv[iEF(t)]dt 
1 Jo 



(11) 



(12) 



and Fit) is defined in Eq[8] Clearly the quasi-energies can 
only become degenerate when they are both equal to zero, and 
we can note from Eq^|that this corresponds, as expected, to 
the destruction of the effective tunneling. 



HI. RESULTS 

To obtain the Floquet quasi-energies for comparison with 
the prediction of Eq[21 the numerical technique described in 
Ref»ii was used. This involves evaluating the unitary evolu- 
tion operator for one period of the field t/(T, 0) and obtain- 
ing its eigenvalues, which are related to the quasi-energies 
via Aj = exp[— iejT], Using this method to obtain the quasi- 
energies, a standard bisection algorithm could then be used to 
find the location of the quasi-energy crossings to a high degree 
of accuracy. 

The dynamical behavior of the system was also examined 
directly by integrating it over long time-periods, with the par- 
ticle initially located in the left quantum dot. To quantify to 
what extent the tunneling between the left and right quantum 
dots was destroyed, the probability that the particle was in the 
left quantum dot (Pl(0) was measured throughout the time 
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FIG. 1 : (a) Quasi-energies for a sinusoidal driving field, of frequency 
a> = 8. Circles indicate exact results, lines the perturbative result 
±(A/2)Jo(E/a>). (b) Localization in the driven system. Spikes in the 
localization are centered on crossings of the quasi-energies. 



evolution. We denote the minimum value of Pi attained dur- 
ing this period to be the "localization", and thus high values 
of localization correspond to the presence of CDT, while low 
values reveal that the particle is able to tunnel from one side 
to the other, and is therefore delocalized. 



A. Sinusoidal driving 

We begin with the most familiar case, when the driving field 
has the form f(f) = cos tot. The procedure outlined in Section 
III B I can be followed straightforwardly, leading to the result 
that: 

i r T 

<<^>r = - exp[iEsm(cot)/cj]dt. (13) 
1 Jo 

By making use of the standard identity: 

DO 

exp[iE sin(iot) / co] = > J,„ (E / co) exp[imiot] (14) 



this expression can be substantially simplified, yielding the fi- 
nal result that e ± = ±(A/2) Jo(E/a>). This reproduces the well- 
known result that for sinusoidal driving CDT occurs when the 
ratio of the field strength to its frequency is equal to a root of 
the Bessel function Jq. In Fig^ the locations of the quasi- 
energies are shown for a fixed frequency co — 8 as a function 
of El oj. It can be seen that the perturbative result works ex- 
tremely well in this regime (high-frequency). Fig^) shows 
the localization produced by the field, as defined above. As 
expected, at the points where the quasi-energies cross the tun- 
neling dynamics of the system is blocked, producing sharp 
spikes in the localization, centered on the crossings. 

To investigate how the crossings move away from these 
points as the driving frequency is reduced, their locations are 
shown as a function of 1 /co in Fig|5^. In accordance with the 
von Neumann- Wigner theorem 815 we can readily see that the 
set of crossings form one-dimensional manifolds. As co tends 
to infinity the crossings occur at the roots of Jq, as predicted 
by the perturbation theory, and this remains a good approxi- 
mation for frequencies as low as co — A. Below this value, 
however, the crossings smoothly drift away from these loca- 
tions, and evolve towards the points A/co - 2n (where n is a 
positive integer), as was seen earlier in Ref. 8 . This limiting 
behavior in the low-frequency regime was also predicted in 
Refpi^, where a similar pattern of crossing-drift was observed 
in an investigation of a related model. The form of FigEt 
immediately suggests fitting the manifolds of crossings with 
quadrants of ellipses: 

where y„ is the n-th root of Jo(y). It can be seen in Fig0 
that this simple parameterization fits the results extremely 
well for the first crossing-manifold, and that the difference be- 
tween the exact location of the crossing and the fitting func- 
tion (Efit/tu - E exac ,/co) never exceeds a value of 0.02. The 
degree of deviation becomes larger as the order of the cross- 
ing increases, but nonetheless is only visible in Fig|5^ for the 
fourth and fifth crossing-manifolds. 



In Fig|4]the localization is plotted as co is reduced from a 
high value towards zero, with E set to hold the ratio E/to on 
a crossing-manifold. For each point the system was evolved 
over 200 periods of the driving field to study how effec- 
tively the field could maintain a localized state. For the high- 
frequency regime, co > A, the localization is excellent at all 
the crossings, with less than 0. 1 of the particle density leaking 
across to the right-side dot during the time-evolution. As can 
be expected, the high-order crossings, which occur at higher 
values of E, can maintain better levels of localization than 
the low-order crossings 8 . This difference becomes more pro- 



nounced as the frequency is reduced, and although the local- 
ization in all cases decays smoothly to zero, the localization 
at the higher-order crossings decays much more slowly. For 
frequencies as low as co = 0.4A, however, the inhibiting effect 
of CDT is still evident for all the crossings, indicating that 
even low-frequency fields may serve a useful role in stabiliz- 
ing electron-leakage from quantum dot devices. 
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FIG. 2: Location of crossings of quasi-energies, in each case the crossings fall on one-dimensional manifolds, 
square-wave driving; (c) triangular driving. Dotted lines indicate the empirical fitting function Eq |15l 



(a) sinusoidal driving; (b) 



first manifold 

second manifold 

third manifold 




FIG. 3: Deviation of the first crossing-manifold from the empirical 
fitting function. Solid line indicates sinusoidal driving; dashed line = 
square-wave driving; dotted line = triangular driving. For the square- 
wave the deviation is smaller than 10" 7 . The dot-dashed line gives 
the deviation for the Fourier expansion of the square-wave, truncated 
at two terms. 



B. Square-wave driving 

Square-wave driving has been considered to a lesser extent 
than the sinusoidal case, although it is also an easily realiz- 
able waveform in experiment. Ref. 20 investigated the case of 
a superlattice driven by a square-wave field, and found that 
for suitable choices of parameters CDT would indeed occur, 
while sinusoidal driving of this system could only produce 
partial CDT 21 . Recently in Ref. 22 it has been shown that in 
a superlattice CDT can only be produced if the crossings of 
the quasi-energies are equally spaced, which clearly does not 
occur for sinusoidal driving. For this reason it is of interest 
to derive the behavior of the quasi-energies for square-wave 
driving to see explicitly how this condition is fulfilled. 

We consider the square-wave driving field, /(f) = 0(f) - 
2@(f - T /2), defined over the interval < f < T. The integra- 
tions required to obtain the quasi-energies may again be done 



FIG. 4: Localization produced by a sinusoidal field, for (E, a>) coor- 
dinates lying on the first three crossing-manifolds. 



straightforwardly, giving the result that: 

A sin (ttE/Ico) 

e± — ± . 

2 tiE/Icj 



(16) 



From this it is immediately clear that the crossings are equally 
spaced as required, being given by the condition E/co - 2n 
where n is a positive integer. In Fig|5^ the quasi-energies ob- 
tained for a frequency of w = 8 are shown in comparison with 
the above result, and it can be clearly seen that the agreement 
is excellent. Below this figure is plotted the localization pro- 
duced by the field, and as for the case of sinusoidal driving, 
the crossings of the quasi-energies correspond to sharp spikes 
in the localization, verifying that CDT is indeed occurring. 

In Fig|2j5 the drifting of the crossings as the frequency is 
reduced is shown. The behavior is strikingly similar to that 
observed for sinusoidal driving, and accordingly we use the 
same functional form d!5l > to fit the crossing-manifolds, with 
the y-intersections now given by y„ = 2n. The fit is so good 
that on this plot no differences can be seen between the exact 
results and the fits. This is corroborated by Fig[3J where the 
deviation from the exact result for the lowest manifold can be 
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FIG. 5: (a) Quasi-energies for a square- wave driving field, of fre- 
quency a> = 8. Circles indicate exact results, lines the perturbative 
result Eq |16l (b) Localization produced by the driving field. 



FIG. 6: (a) Quasi-energies for a triangular driving field, of frequency 
a> = 8. Circles indicate exact results, lines the perturbative result 
Ea ll 81 (b) Localization produced by the driving field. 



seen to be negligible in comparison with the sinusoidal case, 
and within the accuracy of the numerical procedures the fit 
is identical with the exact result. We therefore conjecture that 
this fitting is, in fact, exact for the case of square-wave driving. 
We also show on this plot the result obtained for a bandwidth- 
limited square-wave, obtained by truncating its Fourier expan- 
sion at two terms. We see that the addition of just the second 
term to the sinusoidal driving already reduces the deviation of 
the fit from the exact result considerably. Truncating the se- 
ries at higher points produces steady improvements in the fit, 
strongly supporting the conjecture that the fit is exact when all 
terms are included. 



C. Triangular driving 

We now consider another easily obtainable form of driving, 
the triangular waveform: 



fit) 



1 - At IT 
-3+ At IT 



forO < t < T/2 
T/2<t< T. 



(17) 



For this case a closed form solution can again be obtained for 
the behavior of the quasi-energies, involving the Fresnel sine 
and cosine functions, Six) and C(x). The full expression for 
the quasi-energies is given by: 



[cos(jctt/4)C( V*/2) + sin(x7r/4)5( yfx/2)] (18) 



^Tx 



where x = E/oj. In Fig|6] it can be seen that this func- 
tion is indeed an excellent approximation to the true quasi- 
energies, and that CDT again occurs at the points of quasi- 
energy crossings. The roots of Eq^]may be found numer- 
ically, yielding the result that the first three crossings oc- 
cur when E/co = 2.92519, 7.02525 and 10.9864. Observ- 
ing the behavior of the Fresnel functions 23 reveals that for 
x > 1 they both make small amplitude, decaying oscillations 



about a value of 0.5, which allows the condition for cross- 
ings to be written in the simpler, though approximate, form 
tan(x7r/4) -1. The crossing condition therefore reduces to 
the simple result E/a> - An - 1, as may be seen from the exact 
values given above, which becomes increasingly accurate for 
larger values of E/a>. 

In Fig|2t it can be seen that the crossing-manifolds for this 
form of driving have a similar elliptical form to the previous 
cases. Using the same fitting function (I15> as before, with the 
y-intercepts given by the roots of Eq[D^] gives an accurate de- 
scription of their behavior, as may be seen in Fig [2 Although 
the fit is not as good as for the sinusoidal case, the maximum 
deviation is still less than 0.04. As seen previously, the fit is 
best for the lowest-order manifolds, with small deviations be- 
ing visible in the higher-order manifolds. Nonetheless, in all 
cases the fitting function gives an impressively accurate ap- 
proximation to the true result. 



IV. CONCLUSIONS 

In summary, it has been shown how changing the waveform 
of a periodic driving-field can be used to modify the location 
of the quasi-energies of a two-level system. A procedure has 
been given which relates the waveform explicitly to the quasi- 
energy spectrum, allowing the positions of the quasi-energy 
crossings to be located exactly in the limit of high frequency. 
For various driving fields, including the cases we consider 
here, an analytic form can be obtained for the quasi-energies, 
and in other cases they may be obtained numerically with little 
difficulty. This gives the prospect of designing the waveform 
to create a desired behavior of the quasi-energy spectrum in a 
direct and straightforward way. 

It has also been shown how the positions of the quasi- 
energy crossings drift as the frequency is reduced from the 
high-frequency limit. For the driving fields considered here, 
the crossings fall approximately onto elliptical manifolds, and 
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for the case of square-wave driving it appears that this descrip- 
tion is exact. We have examined this behavior for many other 
waveforms, and we conclude that this form of the crossing- 
manifolds is very general. Using the perturbation theory to 
find the crossings in the high-frequency limit, and then mak- 
ing use of this drifting behavior, allows the positions of the 
quasi-energy crossing to be accurately located in all regimes 
of driving. This gives more flexibility in experiment, as the 
high-field regime may either be difficult to attain, or may in- 
duce undesirable transitions to higher energy levels, breaking 
the two-level approximation. Although the degree of localiza- 
tion that the field can maintain is reduced in the low-frequency 
regime, it can still produce a useful reduction of the leakage 
from quantum dot devices, and thereby enhance their deco- 



herence time, which has many possible applications to the co- 
herent control of mesoscopic systems. 
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